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(57) ABSTRACT 

A method for generating a visual acuity metric, based on 
wavefront aberrations (WFAs), associated with a test subject 
and representing classes of imperfections, such as defocus, 
astigmatism, coma and spherical aberrations, of the subject’s 
visual system. The metric allows choices of different image 
template, can predict acuity for different target probabilities, 
can incorporate different and possibly subject-specific neural 
transfer functions, can predict acuity for different subject 
templates, and incorporates a model of the optotype identifi- 
cation task. 
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PREDICTION OF VISUAL ACUITY FROM 
WAVEFRONT ABERRATIONS 

ORIGIN OF THE INVENTION 

5 

This invention was made by one or more employees of the 
U.S. government. The U.S. government has the right to make, 
use and/or sell the invention described herein without pay- 
ment of compensation, including but not limited to payment 
of royalties. 10 

FIELD OF THE INVENTION 

This invention relates to visual acuity of a human or other 
animal, based on wavefront aberrations associated with the 
animal’s visual imaging system. 

BACKGROUND OF THE INVENTION 

20 

It is now possible to routinely measure the monochromatic 
aberrations of the human eye. However, one cannot yet mea- 
sure the visual acuity that will result from a given set of 
wavefront aberrations. One reason to seek a prediction of 
acuity from aberrations is the possibility of automated objec- 25 
tive measurement of visual acuity, and of automated prescrip- 
tion of sphero-cylindrical corrections. However, it has been 
shown that correcting the spherical and cylindrical compo- 
nents of the aberrations (equivalent to minimizing the RMS 
error of the wavefront) does not provide best acuity. Thus to 
these automated procedures must await a more sophisticated 
metric that can predict acuity from an arbitrary set of aberra- 
tions. 

In the last decade there has been a revolution in measure- 
ment and treatment of visual optical defects. This revolution 35 
has included the development of aberrometers simple enough 
to be used in the clinic, refinement of methods of laser surgery 
for optical correction, and development of various optical 
implants, notably intra-ocular lenses (IOL). In all of these, 4f) 
measurement and interpretation of wavefront aberrations 
(WFAs) has played an important role. They are a simple and 
comprehensive way of describing the state of the optical 
system. In spite of this, there is at present no accepted, reliable 
way of converting WFAs to visual acuity, which is a standard 45 
measure of quality of vision. The WFA Metric allows calcu- 
lation of visual acuity from wavefront aberrations. 

What is needed is an approach, including one or more 
metrics, that allows a prediction of visual acuity, for a human 
or other animal, based on estimated wavefront aberrations 50 
(WFAs) measured or otherwise determined for the test sub- 
ject. Preferably, the approach should allow acuity predictions 
for different optotypes, such as Sloan letters, Snellen e’s, 
Landolt C’s, Lea symbols, Chinese or Japanese characters 
and others. Preferably, the approach should pennit incorpo- 55 
ration of different, possibly subject-specific, neural transfer 
functions. 

SUMMARY OF THE INVENTION 

60 

Uiese needs are met by the invention, which develops and 
applies an optical-based and neural-based metric that allows 
prediction of visual acuity of the subject. For a given choice of 
an optotype set (e.g., Sloan letters), an optical transfer func- 
tion OTF(x,y) is generated, using Zemike polynomials and 65 
the associated Zemike coefficients and a specification of a 
pupil aperture image PA(x,y) for two dimensional coordi- 


2 

nates (x,y) for the subject. A generalized pupil image and 
associated point spread function PSF(x,y) is computed, from 
which an OTF is computed. 

A neural transfer function NTF (x,y) is specified, and a total 
transfer function TTF(x,y) is computed as a product of the 
OTF and the NTF. A proportion correct function P(k) is 
estimated from the neural images and a noise value, using one 
of three or more methods for such estimation. A probability 
criterion P( target) for measurement of visual acuity is speci- 
fied, normally between 0.5 and 0.8. A numerical procedure 
returns a final index value j (final), which is converted to an 
estimate of acuity using a standard logMAR calculation. The 
output of the logMAR computation is a WFA metric that 
provides an estimate of visual acuity for the subject. 

The metric(s) developed here is designed to predict symbol 
acuity from wavefront aberrations. One embodiment of the 
metric relies on Monte Carlo simulations of a decision pro- 
cess and relies on an ideal observer, limited by optics, neural 
filtering, and neural noise. A second metric is a deterministic 
calculation involving optics, symbols, and a hypothetical 
neural contrast sensitivity function CSF. 

A WFA Metric is an algorithm for estimating the visual 
acuity of an individual with a particular set of visual wave- 
front aberrations (WFAs). The WFAs represent arbitrary 
imperfections in an optical system, and can include low order 
aberrations, such as defocus and astigmatism, as well as high 
order aberrations, such as coma and spherical aberration. 
WFAs can now be measured routinely with an instrument 
called an aberrometer. In modern practice, the WFAs are 
represented as a sum of Zemike polynomials Z(x,y), each 
multiplied by a Zernike coefficient. A typical measurement 
on the eye of a subject will consist of a list of about 16 
numbers, which are the coefficients of the polynomials. The 
WFA Metric converts the list of numbers into an estimate of 
the visual acuity of the subject. If changes are planned to the 
WFA of the subject (through surgery or optical aids) the 
predicted change in visual acuity can be calculated. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates an embodiment of a structure for compu- 
tation of a WFA metric according to the invention. 

FIG. 2 illustrates the ten Sloan letters, expressed in a sans 
serif font, that provide one of the optotype sets that can be 
used with the invention. 

FIGS. 3A-3E illustrate steps in creation of the OTF. 

FIG. 4 graphically illustrates a representative radial neural 
transfer function. 

FIG. 5 is an embodiment of a procedure for evaluation of a 
probability correct index for one size optotype. 

DESCRIPTION OF THE INVENTION 
0. Notation and Terminology 

In this presentation, an “image” refers to a finite discrete 
digital image represented by a two-dimensional array of inte- 
gers or real numbers. It has a width and height measured in 
pixels. Where the size is specified it will be given as a list 
{rows, columns}. The image has a resolution measured in 
pixels/degree. The pixel indices of the image are x (columns) 
and y (rows). Images will usually be an even number of pixels 
wide and tall. If the image size is {2h,2h}, then the indices x 
and y each follow the sequence { -h, . . . , 0, . . . , -h- 1 } . This 
places the origin of the image at the center. An image may be 
written with explicit row and column arguments A(x,y), or 
without the coordinates as A. 
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In this presentation, a dft refers to a two-dimensional finite 
discrete digital array of complex numbers representing a Dis- 
crete Fourier Transform (DFT). It has a width and height 
measured in pixels. Where the size is specified it will be given 
as a list {rows, columns}. A dft has a resolution measured in 
pixels/cycle/degree. The pixel indices of the dft are a (col- 
umns) and v (rows). Dfts will usually be an even number of 
pixels wide and tall. If the dft size is {2h,2h}, then the indices 
u and v each follow the sequence {0, . . . , h-1, -h . . . , — 1 }. 
This places the origin of the dft at the first pixel. This is the 
conventional ordering of indices in the output of the Fast 
Fourier Transform (FFT) operator. The FFT is a particular 
algorithm for implementation of the DFT. In the body of this 
document we refer to the DFT, but this will usually be imple- 
mented by the FFT. 

In this presentation, vectors will be written with one sub- 
script Aj., and matrices will be written with two subscripts 
A Jk , where the first subscript indicates the matrix row. Fre- 
quently, we will deal with vectors or matrices whose elements 
are images, in which case the image coordinates x,y are 
omitted. 

We make use of the notation A:B to indicate Frobenius inner 
product of two matrices 

A(x, y)B(x, y ) 


This is useful to describe a sum over pixels of the product of 
two images. The modulus or norm of an image is given by 

||A|[WAiA 

1 . Inputs and Output 

The WFA metric has four inputs. A first input is a set of 
wavefront aberrations, represented as a weighted sum of 
Zemike coefficients z„(x,y). A second input is a set of opto- 
types, represented in a standard graphic format, such as a font 
description, a set of raster images, or graphic language 
descriptors. One example set of optotypes is the Sloan font for 
the letters {C, D, FI, K, N, O, R, S, V, Z}, a set often used in the 
measurement of acuity). A third input is a set of templates, 
equal in number to the number of optotypes in the set. By 
default, the templates are derived from the optotype set and 
are not a distinct input. A fourth input is a set of parameters, 
some of which may have default values that are permanently 
stored within the program. Some parameters may be changed 
on every calculation of the metric, while others are unlikely to 
be changed often. The parameters are described throughout 
this description. 

A single output, the visual acuity, is expressed as a decimal 
acuity or log of decimal acuity (logMAR). An overall system 
structure is shown in FIG. 1. 

2. Overview of the Algorithm 

a. Generate the Optical Transfer Function (OTF) 

b. Generate the Neural Transfer Function (NTF) 

c. Generate the Total Transfer function (TTF) 

d. Define the Proportion Correct function P(size) 

e. Find the size for which P(size)«P, OJg( , r 

Each of these steps is described in detail in the following. 

3. Select a Set of Optotypes 

The optotypes are a set of graphic symbols that the human 
observer is asked to identify in the course of an acuity test. 


4 

Examples are Sloan letters, Snellen e’s, Landolt Cs, Lea 
symbols, Chinese or Japanese characters, or other pictograms 
of various sorts. Each optotype set will have a fixed number K 
of elements, and a defined size specification. 

5 By way of example, the optotype set used here is the Sloan 
letters {C, D, H, K, N, O, R, S, V, Z}, with K=1 0. These letters 
are shown in FIG. 2. Each Sloan letter lias a stroke width 
MAR, expressed in minutes of arc of visual angle, and each 
letter is 5 MAR tall by 5 MAR wide. The size specification 
10 used here is Log 10 MAR, expressed as 

log MAR(mar)=log l0 (mar) 

4. Determine the Usable Range of Optotype Sizes 
15 

The usable range will be limited by the resolution and size 
of the PSF image. As discussed below, these are determined 
by the pupil size, the wavelength (X), and the pupil magnifi- 
cation (m). If the PSF image has a width of r, expressed in 
20 pixels, and d in degrees, the smallest stroke-width possible is 
one pixel, or 

7 60^ 

25 logM/W? mi „ =log 10 ^ — J 

The largest stroke-width will be one fifth width of the largest 
character, which will be one half the width of the PSF image; 
3Q a margin is required to accommodate blur and to avoid wrap- 
around so that 

iog.i«i; m „=iog 10 (6rf) 

It is sometimes convenient to adopt a positive integer index 
35 that corresponds to size. One example is computing logMAR 
in steps of 1/20. In that scheme, the minimum and maximum 
indices would be 

index„,,„=Ceiling(20 log MAR min ) 

40 

index„ jax =Floor(20 log MAR max ) 

The index 1 then extends from 1 to l max =index max -index m! „+ 
1 , and log MAR is given by 

45 

/ + index min - 1 

log/tMR = — , where 1=1,... , l max 

5Q Using the default parameters, the PSF image will have a width 
of 256 pixels, and a width of 0.81 5525 deg. With these values 

i ndcx rnirl - 1 4 
55 ui dex ™ m = 13 

The size index 1 will have values between 1 and l max =28 for 
this example. 

5. Generate the Optical Transfer Function (OTF) 

60 

The mathematical operations required to generate an opti- 
cal transfer function (OTF) from a set of Zernike polynomials 
are well known. Graphs of the results at several stages are 
shown in FIGS. 3A-3E. 

65 a. Create the Pupil Aperture Image PA(x,y). The image is of 

size {2h,2h}, where h is the half width of the pupil 
aperture image, expressed in pixels, 
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PA(x , v) = I if y/ x 2 + y 2 j h < 1, x and y integers e{-/i,...0, 
0 otherwise 


■ h- 1 } 


6. Generate the Neural Transfer Function (NTF) 

a. The Radial Neural Transfer Function RNTF(u,v) is a 
two-dimensional real dft given by 


b. From the set of Zemike coefficients C={c 0 , c x , c 2 , . . . , 
c^} (expressed in microns), create a discrete digital 
image of the Wavefront Aberration Image WA(x,y), with 
image size {2h,2h}. If the Zemike polynomials z„(x,y) 
are identified by single index (the mode) n=0, . . . N; and 
if the c k are the coefficients of the individual polynomi- 
als, then 


RNTF(u, v) = gam|exp -|y- j - lossexp ) | 


where gain, f 0 , f 1; b, and loss are parameters. An example of 
this function is shown graphically in FIG. 4. 

b, The Oblique Effect Filter OEF (u,v) is a two-dimensional 
real dft given by 


15 


WA(x, y) = z c„Zn(x, y) 


We make use of the standard form of the Zernike polynomials 
as defined by Thibos, 2002, Jour. Of Optical Society of 
America. 

c. Compute the Generalized pupil image GP(x,y) 


20 


{ ( f - corner!! 

OEF{u , v) = OEF{f , hi = I - 1 1 - cxp| -j JJsin(2(?) if / > comer 


= 1 otherwise 


/UlF+v 2 


GP(x, y) = PA(x , y)epx WA(xv,) 


25 


6=arctan(«,v) 


where X is the wavelength of light in rnn used to illuminate the 
optotype set. 30 

d. Pad the image on the left and top with zeros to create tut 
image of size {2hm, 2hm}. The parameter m is the pupil 
magnification. 

e. Compute the Point Spread Function PSF(x.y) 

, 35 

PSF(^)=IDFT[GP(x,y)]l 2 

where DFT is the Discrete Fourier Transform operator. 

f. Normalize the PSF. 


40 


PSF(x, y) ■■ 


PSF(x, y) 

II PSF(x, v)|| 


g. Compute the Optical Transfer Function OTF(u.v) 

OTF(u,v)=2AmDFT[PSF(;t,j)] 

This result is a complex image of size {2hm, 2hm}. 

The height and width of the PSF image in degrees of visual 
angle is given by 


45 


50 


d = - 


/i360A10" 6 
px 


where p is the pupil diameter in mm. The height and width of 
the PSF image in pixels is given by 


where corner and slope are parameters. 

c. Compute the Neural Transfer Function NTF(u,v), a two- 
dimensional real dft given by 

NTF( a, v)=RNTF( u, v)OEF(u, v) 

7. Generate the Total Transfer Function (TTF) 

The Total Transfer Function is given by 

TTF(u,v)=OTF(a, v)NTF(«, v) 

8. Define the Proportion Correct Function P(k) 

The steps in evaluation of the P(k) function are as follows, 
and are diagrammed in FIG. 5. 

a. Given a size index 1 , create K optotype images O fc (x,y). 
This may be done by rendering images from a graphic 
description, or the images may be pre-computed. Each 
image is of size {r,r}. See above for a definition of the 
optotype size index 1 . 

b. Create the K Neural Images S fc (x,y) by computing the 
DFT of the each optotype image O^, multiplying by the 
TTF, and taking the inverse DFT, 

S*=IDFT[DFT[O t ]TTF] 

where DFT is the DFT operation and IDFT is the Inverse 
DFT operation. 

c. Create the K template images T k . By default, these are 
identical to the Neural Images S k . 

d. Compute the normalized templates. Each template is 
divided by its norm, equal to the square root of the sum 
of the squares of all its pixels. 


r=2hm 


where h is a half-width. The resolution of the PSF image in 60 
pixels/degree is 


IIA II 


r 2 nmp 

V ~d~ 360AKT 6 


e. Compute the matrix of normalized template cross-cor- 


relations W. 


jjc 
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f. Create an array of cross-correlations between each neural 
images and each template. Note that the row indexes the 
neural image and the column, the template. 

%,t = S,:T fe 5 

g. At this point two or more methods are available, which 
we identify as methods 1 and 2. 

Method 1 . 

i. Subtract each value from the main diagonal entry in the 
same row, and divide by a factor that includes the param- 1 0 
eter cr (default value** 1 ). There are two possible versions 
of a matrix D, identified by subscripts 1 and 2. 


Dl.iL = - 


r jj ~ R iJj 

"V 1 - w h 


15 


ii. The probability correct for optotype j is given by 


20 


R i.j = rv«n F(t-D UJJL )dt 

J -°° k*j 


25 


where f(t) and F (t) are probability density function and cumu- 
lative probability 
Method 2 


Dl,j,k '■ 


R jJ - r M 
crsfl^ 1-Wj, k 


P2j = n 

k*j 


30 


35 


The final value of P is given by 


40 


1 K 

p= p k 


45 

9. Find the Size for which P m P, arget 

The parameter P te , is the criterion probability for mea- 
surement of visual acuity. It is usually set to a value between 
0.5 and 0.8. This value will depend upon the number K of 50 
optotypes and must be greater than 1/K (the probability of 
getting the right answer by guessing). For the Sloan letters, a 
default value P Mfger =0.55 is used. Various efficient iterative 
procedures may be used to locate the value of size for which 
PM 3 target- Here we describe the method of bisection, though 55 
other methods may be used. 

how = l 


Plow PQlow) 

Phigh = P^high) 


If Ii..,,!, -liiw — 1, exit and return / 


, * high ‘-tow . . 

'-final = •-low ^ \Pt — Plow ) 

Phigh - Plow 


Lid = Round 


Lgh "P llow 1 


P mid PQmid) 


If Pw<P 


target ’ 


1 


=1 




dlow) 


otherwise 

•high •mid 


Phigh PQhigh) 

Go to begin loop 

The returned value of \ flnal can then be converted to an acuity 
in logMAR using the Equation above. This is the output of the 
WFA Metric. 

FIG. 5 illustrates a sequence of steps of a procedure for 
practicing the invention. In step 51, an OTF is generated. In 
step 52, an NTF is generated and is multiplied by the OTF, to 
form a TTF (step 53). In step 54, a set of optotypes is and is 
subjected to a DFT process, in step 55. In step 56, the pro- 
cessed optotypes are used to form images S, of the optotypes. 
In step 57, the images S y are used to create a set of templates 
T, t , and normalized templates T** are created in step 58. 
Cross-correlations R ; i of the images S ; and the normalized 
templates T^* are formed, in step 59. In step 60, cross-corre- 
lations W k of the normalized templates T^*. Normalized dif- 
ference matrices D yjt are formed from the cross-correlation 
matrix R ; A . are formed in step 61, using information from the 
cross-correlations W k . and a statistical parameter a, in step 
62. In step 63, a probability P associated with measurement of 
visual acuity is computed. 

10. Unique Features of the WFA Metric 

The WFA metric is the only known metric to compute 
acuity from wavefronts that: 

(i) incorporates a model of the optotype identification task 

(ii) can predict acuity for different target probabilities 

(iii) can predict acuity for different optotypes 

(iv) allows user specification of optotypes 

(v) can incorporate different and possibly subject-specific 
neural transfer functions 

(vi) can predict acuity for different subject templates 
The template matching algorithm that is fundamental to this 
metric may have other uses in predicting performance in 
identification tasks. 


Default 

Parameter value Unit Definition 


K 

k 

L 

556 

nm 

number of optotypes 
index of optotpye, 1 , . . . , ! 
wavelength 

P 

5 

mm 

diameter of pupil 

h 

64 

pixels 

half width of pupil image 

m 

2 


magnification 

d 

derived 

degrees 

size of the PSF image 


begin loop 
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-continued 



Default 



Parameter 

value 

Unit 

Definition 

r 

2 m h 

pixels 

size of the PSF image 

V 

r/d 

pixels/deg 

resolution of PSF image 

comer 

13.5715 

cycles/deg 

oblique effect parameter 

slope 

3.481 


oblique effect parameter 

gain 

3.149614 


NTF parameter 

loss 

0.9260249 


NTF parameter 

fo 

35.869213 


NTF parameter 

fi 

5.412887 


NTF parameter 

b 

1 

1.064181 


NTF parameter 
optotype size index 

WA 


image 

wavefront aberration 

PSF 


image 

point spread fimction 

PA 


image 

pupil aperture 

GP 


complex image 

generalized pupil 

TTF 


dft 

total transfer fimction 

NTF 


dft 

neural transfer fimction 

OTF 


dft 

optical transfer fimction 

OEF 


dft 

oblique effect transfer fimction 

O* 



optotype with index k 

S* 



neural image of optotype 

T* 



template with index k 

T* 



normalized template with index k 

w,.* 



cross-correlation between 

^j,k 



normalized templates 
cross-correlation between 



normalized templates 
and neural images 





template response distribution 




means 

P/ 


probability 

probability correct for optotype 
with index k 

P 


probability 

probability correct for optotypes 
of one size 

a 



noise standard deviation 

Y targl , t 

0.55 


criterion proportion correct 

mar 


minutes 

optotype stroke size 

x,y 



image pixel coordinates 

u,v 



dft pixel coordinates 

j,k 



row, column indices of matrices 

C n 



coefficient of Zernike 
polynomial n 

Z„ 



Zemike polynomial n 


What is claimed is: 

1. A method of predicting visual acuity of an imaging 
system, the method comprising: 

providing a set of K optotypes, numbered k=l,. . . , K 
(KS2) of a specified size, to be used to establish a 
reference set of symbols to estimate visual acuity, and 
providing a description or image Tj.(x,y) of each opto- 
type, dependent upon location coordinates (x,y) associ- 
ated with the optptype image produced by the imaging 
system; 

constructing a Wavefront Aberration Image WAI(x.y) that 
manifests selected image aberrations associated with the 
imaging system; 

computing a generalized pupil image GPI(x,y), with non- 
zero values confined to within a pupil aperture associ- 
ated with the imaging system, that manifests the selected 
image aberrations and that is dependent upon at least one 
wavelength X of light with which the image is viewed; 
computing a Point Spread Fimction PSF(x,y) that is an 
absolute value squared of a Discrete Fourier Transform 
(DFT) of the image GPI(x,y); 
computing a Normalized Point Spread Function NPSF(x, 
y), proportional to the Point Spread Function PSF (x,y), 
whose norm is 1: 

computing an Optical Transfer Function OTF(u,v), 
expressed as a fimction of spatial frequency indice (u,v) 
in a Fourier transform plane, that is a Discrete Fourier 
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Transform of the Normalized Point Spread Function 
NPSF (x,y), multiplied by a value h of original image 
size and multiplied by an image magnification index m; 
generating a Radial Neural Transfer Fimction RNTF(u,v), 
5 where RNTF(u,v) satisfies the following conditions: (i) 
RNTF(u,v) is a continuous, non-negative fimction of a 
spatial frequency variable f=[u 2 +v 2 ] 1/2 ; (ii) RNTF(u,v) 
has at least one value, f=f(max), for which NRTF(u,v) is 
a maximum; (iii) for 0<f<f(max), RNTF(u,v) is mono- 
10 tonically increasing in f; (iv) for f>f(max), and RNTF 
(u,v) is monotonically decreasing toward 0; 
computing a Neural Transfer Function NTF (u,v) as a prod- 
uct of the Radial Neural Transfer Function RNTF (u,v) 
and an Oblique Effect Filter fimction OER(u,v) that 
15 compensates for viewing angle of the-original optotype 
image; 

generating a Total Transfer Function TTF(u,v), defined as 
a product of NTF (u,v) and the Optical Transfer Function 
OTF(u,v); 

20 generating a Proportion Correct index PC(j), which pre- 
sents a probability associated with a correct optotype 
that would be identififed by a subject having the Wave- 
front Aberration Image WAI(x,y) for at least one opto- 
type(j). 

25 2. The method of claim 1, wherein said step of constructing 

said Wavefonn Aberration Image WAI(x,y) comprises: 
providing a Pupil Aperture Imaging characteristic function 
PA(x,y) that has a first value substantially equal to 1 
within a pupil aperture of a test subject and has a second 
30 value substantially equal to 0 outside the pupil aperture 
of the test subject, where (x,y) are location coordinates; 
providing a set of Zemike polynomials Z„(x,y) and asso- 
ciated coefficients c n , and creating a discrete digital 
image a Wavefront Aberration Image WAI, defined as an 
35 error sum 

WAI(x,yy=2„c„Z n (x,y); 

where the coefficients c„ are chosen to minimize a computed 
error between an image provided by said imaging system and 
40 the error sum; 

computing a generalized pupil image GP, defined as 

G?(x,y)=PA(x,y)exp{2iWAl(x,y)/QJ1000)}, 

45 where X refers to a wavelength used to illuminate the opto- 
types; and 

padding said image of each of said optotype on at least two 
adjacent edges of said image with zeroes to create an 
image of size (Ax,Ay)=(mh,mh), where h represents an 
_ f) original image size, measured in a selected direction, 
and m is a multiplier index. 

3. The method of claim 1, wherein said Oblique Effect 
Filter OEF(u,v) is determined by a procedure comprising 
forming 

55 

OEF(u, v) = OEFl(f, 9) = 1 - {1 - exp{((corner) - /)/ (slope )}}sin20 
if / > (corner) 

=1 otherwise 


f=(u 2 +v 2 ) 1/2 

0=arctan(«,v), 

where comer and slope are parameter. 
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4. The method of claim 1, wherein said process of gener- 
ating said Proportion Correct index PC(j) comprises: 

computing a neural image S^(x,y) (j=l , . . . , K) as an inverse 
DFT of a product of said Total Transfer Function TTF 
(u,v) multiplied by a DFT of an optotype image 0 ; (x,y); 
computing a cross-correlation matrix between each of the 
neural images S . and each of said templates T t , defined 
as 

R, J rS/ x ,y):T i (x,y); 

computing a cross-correlation matrix of each normalized 
template with each normalized template, W / (t = T^T*.; 
computing and normalizing a difference of matrix entries, 
Rj -R ,ij„ between each diagonal entry R, and each entry 
Rj k in a corresponding row of the matrix {R ; j.}, to fonn 
a normalized difference matrix D ^.; 

D,^={R, v -R,^}/{a{l-W, Jt 2 } 1 ' 2 }, 
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where a is a statistical value that is provided or computed; 
and 

estimating a probability for estimation of a correct opto- 
type by a subject by computing a probability value that is 
5 either (i) a product of functions of the matrices D or (ii) 

an integral, with an integrand equal to a normal statisti- 
cal function, multiplied by the product of the functions 
of the matrices T> Jk . 

5. The method of claim 1, wherein said Radial Neural 
to Transfer Function RNTF(u,v) is defined as 

RNTF(u,v)=(gam){exp[-(f/f 0 ) i ’-(loss) exp[-(f f L ) 2 ]}, 

where f=(u 2 +v 2 ) 1/2 is a spatial frequency and (gain), (loss), 
b, f 0 and f : are selected parameter values. 

15 6. The method of claim 1, wherein said set of optotypes 

includes at least one of: 

Sloan letters; Snellen E’s; Landolt C’s; and Lea symbols. 



